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Abstract. We present a study of Ar abundances in 15 damped Ly a systems (DLAs) in the redshift interval 
2.3 < z < 3.4. The sample includes 4 UVES/VLT measurements of An column densities presented here for 
the first time, together with 6 measurements and 5 upper/lower limits collected from the literature (UVES/VLT 
and HIRES/Keck observations). The majority of DLAs show significant underabundances of Ar relative to other 
a-capture elements with common nucleosynthetic origin. We show that neither dust depletion nor intervening 
Hn regions inside DLAs offer a viable justification to these underabundances. A natural explanation is found 
in the framework of photoionization models of Hi regions embedded in an ionizing continuum with varying 
spectral distribution. At z ~ 2.5 the observed Ar deficiencies are large, [Ax/a] ~ —0.6/ — 0.8 dex, suggestive of 
a hard, QSO-dominated spectrum. At z > 3 the deficiencies are instead small, suggestive of a soft, stellar-type 
spectrum, though more data are needed to generalize this high- z result. Should the change of Ar abundances with 
redshift be governed by the evolution of the UV stellar emission internal to DLAs, a synchronization of the star 
formation in DLAs would be required, with a strong stellar emission at z > 3, but weak at z < 3. This variation 
seems inconsistent with the weak signal of evolution indicated by abundance studies of DLAs. More likely, the 
change of Ar abundances is induced by the evolution of the UV metagalactic continuum, in which case the UV 
emission internal to DLAs must be small (i.e. DLAs should have modest star formation rates) and the external 
background must become softer at z > 3. The former requirement is consistent with the modest evolution of DLAs 
abundances and the lack of Lya and Ha emissions associated with DLAs. The latter requirement is consistent 
with the observed evolution of Si iv/C IV ratios in the IGM, the claims of high escape fraction of UV photons from 
Ly-break galaxies at z > 3, and the recent finding that the Hell re-ionization seems to occur between z ~ 3.4 and 
z ~ 3. Comparison with results from local interstellar studies indicates that Ar abundances can be used to trace 
the evolution of the ionization history of the universe down to z = 0, where [Ar/a] ~ —0.2 dex. We predict a rise 
of Ar abundances in the redshift range from z ~ 2.3 to z = 0, at the epoch at which the metagalactic field of 
galaxies overcomes that of quasars. 
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1. Introduction 

The abundance of argon in diffuse gas can be measured 
from absorption line spectroscopy of the An 104.8 and 
106.6 nm resonance transitions. Early Copernicus obser- 
vations indicated that argon is moderately depleted in the 
local interstellar gas (Meyer 1989 and refs. therein), a re- 
sult confirmed by recent IMAPS (Sofia & Jenkins 1998; 
hereafter SJ98) and FUSE (Jenkins et al. 2000; Lehner et 
al. 2002) observations. Interstellar depletions are gener- 
ally attributed to the incorporation of atoms into dust 
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grains and tend to be correlated with the gas column 
density (Jenkins 1987; Savage & Sembach 1996). It is 
hard to establish whether such a correlation, suggestive 
of dust, exists or not for Ar because An lines saturate 
when iV(Hi) > 10 20 atoms cm -2 at the solar metallicity 
characteristic of the local interstellar gas. In fact, SJ98 
provided theoretical arguments to show that Ar depletion 
is unlikely due to dust since Ar atoms have a low proba- 
bility of being incorporated into dust grains. The Ar de- 
ficiency is most likely due to ionization effects since the 
ratio between photoionization and recombination rates is 
typically one order of magnitude larger for Ar I than for 
Hi (Fig. 3 in SJ98). Owing to this property, Lyman con- 
tinuum photons with hv > 15.76 eV, the An ionization 
threshold, will be extremely efficient in ionizing Ar I but 
not Hi, if they are able to leak through a Hi layer. This 
characteristic of Ar can be used to probe ionization con- 
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Table 1. Column densities" and relative abundances of argon in Damped Ly a systems. 
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log iV(Hl) 
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[Ar/O] 


[Ar/S] 


[Ar/Si] 


A b 

A [Ar/Sil 


References 


0528-2505 


2.812 


21.11+0.04 


14.25+0.01 




-0.55+0.02 


-0.63+0.04 


-0.04 


This work; Pcroux et al. (2003) 


0841+129 


2.375 


21.00+0.10 


13.46+0.10 




-0.70+0.10 


-0.74+0.10 


-0.00 


This work; Centurion et al.(2003) 


0841+129 


2.476 


20.79+0.10 


13.34+0.12 




-0.57+0.12 


-0.55+0.14 


-0.01 


This work; Centurion et al.(2003) 


2343+1232 


2.431 


20.18+0.05° 


13.22+0.05° 




-0.78+0.05° 


-0.84+0.05° 


— 


This work; D'Odorico et al.(2002) 


0000-263 


3.390 


21.41+0.08 


14.02+0.03 


-0.03+0.03 


+0.00+0.04 


+0.00+0.04 


-0.00 


Molaro et al. (2000,2001) 


0201+365 


2.462 


20.38+0.05 


14.08+0.03 




-0.53+0.03 


-0.42+0.03 


-0.20 


Prochaska et al. (2002a) 


0347-383 


3.025 


20.63+0.01 


13.99+0.02 


-0.43+0.02 




-0.21+0.02 


—0.00 


Levshakov et al.(2002) 


0741+4741 


3.017 


20.48+0.10 


13.13+0.02 






-0.18+0.02 




Prochaska et al. (2002a) 


1759+7539 


2.626 


20.61+0.01° 


13.68±0.03 c 




-0.71+0.13° 


-0.68+0.04° 




Prochaska et al. (2002b) 


2243-6031 


2.330 


20.24+0.02° 


13.35+0.05° 




-0.59+0.04° 


-0.54+0.04° 


-0.00 


Lopez et al.(2002) 


0001 


3.000 


20.70+0.05 


<13.38 


< -0.16 




< -0.03 




Prochaska et al. (2002a) 


0336-0142 


3.062 


21.20+0.10 


>13.94 




> -0.37 






Prochaska et al. (2002a) 


0930+2858 


3.235 


20.30+0.10 


<12.96 






<+0.11 




Prochaska et al. (2002a) 


1426+6039 


2.827 


20.30+0.04 


<13.43 






< -0.36 




Prochaska ct al. (2002a) 


2344+12 


2.538 


20.36+0.10 


<13.26 






<+0.12 




Prochaska ct al. (2002a) 



a Column densities are expressed in atoms cm -2 . 

Additive dust correction term for [Ar/Si] estimated using the set of parameters Sll(Vladilo 2002) and assuming zero dust depletion for Ar. 
° Sum of velocity components observed in An; the fraction of TV(Sill) measured in such components over total A^(Siir) has been used to scale down 
the total AT(Hl); components 6 and 7 in Prochaska et al. (2002b) have been considered for the DLA in Q1759+7539; components 3, 4 and 11 in 
Lopez et al. (2002) for the DLA in Q2243-6031. 



ditions not only in local interstellar clouds, but also in 
clouds of external galaxies. In particular, with high reso- 
lution spectroscopy of QSO sources we can study clouds in 
intervening galaxies detected as absorption line systems. 
Damped Ly a systems (DLAs), are thought to arise in 
these intervening galaxies (Wolfe et al. 1995) and offer 
the possibility to probe ionization conditions at high red- 
shift from the study of Ar abundances. At z > 2.3 the Ar I 
lines can in fact be detected in the optical range, where the 
collecting power of 10-m class telescopes can be exploited 
to obtain high quality QSO spectra. Thanks to the low 
metallicity typical of DLAs (Pettini et al. 1999), the Ar 
resonance lines can be unsaturated even at the high col- 
umn densities, AT(Hi) > 2 x 10 20 cm~ 2 , typical of these 
systems. Therefore, DLAs can yield information comple- 
mentary to that obtained from studies of local ISM, where 
the An lines become readily saturated. Most important, 
if Ar abundances in DLAs are indicative of ionization con- 
ditions, as they arc in the local ISM, their study may yield 
important constraints on the nature and relative impor- 
tance of ionizing sources at high redshift, a topic of par- 
ticular interest in cosmology for its effects on the physical 
state of the intergalactic medium and for its link to the 
star formation history of the universe. 

The first detection of An in a DLA system was ob- 
tained for the z = 3.39 absorber towards QSO 0000- 
26 observed with UVES/VLT (Molaro et al. 2001). The 
Ar abundance was found to match that of other a- 
capture elements measured in the same system, suggest- 
ing a low degree of ionization of the gas. Recently a few 
other Arl detections in DLAs have been obtained with 
UVES/VLT (Levshakov et al. 2002, Lopez et al. 2002) 
and HIRES/Keck observations (Prochaska et al. 2002a,b). 
From these works there was no evidence for a general un- 
derabundance of Ar in DLAs, with the notable excep- 
tion of the z = 2.62 system towards QSO 1759+75, an 
absorber characterized by complex ionization conditions 
(Prochaska et al. 2002a). 



From the theoretical point of view, the Ar ionization 
fractions in DLA systems have been estimated by several 
authors, assuming different models of radiation field and 
ionization structure of the absorbers (Izotov et al. 2001, 
Vladilo et al. 2001, Prochaska et al. 2002a). At variance 
with other metal species observed in DLAs, the fraction 
of Ar I is quite sensitive to the adopted model. This makes 
An an important discriminant of ionization conditions 
in DLAs. Motivated by this possibility, we have actively 
searched for new Ar I detections in DLAs. In this paper 
we present and discuss the results of this search. The new 
measurements are described in Section 2, while in Section 
3 we review the general properties of Ar abundances in 
DLAs combining these new data with previous measure- 
ments. We show that a significant underabundance of Ar 
is common among DLAs, contrary to previous belief, and 
we discuss this result in the framework of photoioniza- 
tion equilibrium models. In Section 4 we present evidence 
for redshift evolution of the Ar abundances in DLAs and 
we briefly summarize Ar interstellar abundances at z = 0. 
Finally, in Section 5 we discuss these findings in the frame- 
work of current evolution scenarios for DLAs and for the 
metagalactic radiation field. The results are summarized 
in Section 6. 

2. The data 

The observations were obtained with the Ultraviolet 
Visual Echelle Spectrograph (UVES) fed by the Kueyen 
telescope (unit 2 of the ESO VLT). This spectrograph 
is particularly well suited for observing the spectral re- 
gion in the proximity of the ultraviolet atmospheric cutoff 
(Dekker et al. 2000) allowing the Ar I lines to be observed 
at redshifts z > 2.3. The slit widths were set to 1 arcsec, 
yielding a mean resolving power R = AA/A ~ 4.2 x 10 4 . 
The targets, the basic data for the DLAs investigated, and 
the new measurements of An column densities, V(An), 
are listed in the first 4 rows of Table 1 . A full presentation 
of the data reduction is given in a parallel paper focussed 
on a study of nitrogen abundances in DLAs (Centurion 
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Fig. 1. Absorption profiles of An and S n lines (see labels) 
for the DLA system at z = 2.8120 (zero velocity) in QSO 
0528-2505. 
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Fig. 2. Absorption profiles of An and Sn lines (see la- 
bels) for the DLA system at z = 2.3745 (zero velocity) in 
QSO 0841+129. 



et al. 2003). More details on UVES studies of the sys- 
tems at z = 2.812 in QSO 0528-2505 and z = 2.431 in 
QSO 2343+1232 can be found in Peroux et al. (2003) and 
D'Odorico et al. (2002), respectively. The column densities 
have been derived with the Voigt-profile fitting routines 
FITLYMAN (Fontana & Ballester 1995) implemented in 
the MIDAS software package. The oscillator strength for 
the transition An A104.8nm, / = 0.2570, was taken from 
Federman et al. (1992), the same adopted for the local ISM 
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Fig. 3. Absorption profiles of Ar I and S n lines (see labels) 
for the DLA system at z = 2.4762 (zero velocity) in QSO 
0841+129. 
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Fig. 4. Absorption profiles of Ar I and S II lines (see la- 
bels) for the DLA system at z = 2.4312 (zero velocity) in 
QSO 2343+1232. 
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by SJ98 and in previous work on DLAs by Molaro et al. 
(2001) and Levshakov et al. (2002). Oscillator strengths 
for other a-capture elements were taken from Morton 
(1991), except for: Sin 1304, 1526 (Spitzer & Fitzpatrick 
1993) and Sin 1808 (Bergeson k Lawler 1993). 

The absorption profiles of An lines are shown in Figs, 
n 121 121 an d El Care has been taken in assessing possible 
contamination of the Ar I profiles due to Lya interlop- 
ers, which would tend to increase A^(Ari). Comparison 
with the radial velocity profiles of low-ionization species 
falling inside and outside the Lya forest has not revealed 
evidence of contamination. In the figures we show exam- 
ples of the profiles of the Sn and/or Sin lines used in 
the present analysis. The agreement with the radial ve- 
locity profiles of An lines is generally excellent. A slight 
misalignment between the Ar I and S n lines is only ob- 
served in the z = 2.3745 system in Q 0841+129. We have 
carefully checked that this is not due to an error in the 
wavelength calibration. Some weak contamination could 
be responsible for this shift, in which case the derived Ar 
column density could be slightly overestimated, probably 
within the large error bar of this measurement. This would 
not affect the conclusion that Ar is significantly under- 
abundant in this system. Care has been taken in assessing 
systematic errors due to the uncertainty of the continuum, 
which are largest for the two DLAs in Q 0841+129 (Figs. 
Eland Even in these two cases the A^(An) error bars 
are sufficiently large to encompass the continuum errors. 

The 4 new measurements presented here repre- 
sent a significant increase over the 6 A^An) measure- 
ments in DLAs previously obtained with UVES/VLT or 
HIRES/Keck observations, also shown in Table 1. When 
also upper and lower limits of A^(An) arc considered (last 
5 rows of Table 1), the combined sample sums up to a 
total of 15 DLAs. Care has been taken in comparing the 
abundances derived by different authors. The / value of 
the An 104.8 nm line adopted by us is 2% smaller than 
the updated Morton (1991) value adopted by Jenkins et 
al. (2000) and Lehner et al. (2002) for ISM studies and 
by Prochaska et al. (2002a, b) for DLAs studies. The value 
is 5% larger than the old Morton (1991) value adopted 
by Lopez et al. (2002). We did not attempt to correct for 
these differences, given the presence of some saturation in 
the lines. The differences are in all cases well within the 
measurement errors quoted by the authors. 

For the sake of comparison with previous work, in Fig. 
we plot the [Ar/Fe] versus [Fe/H] abundance ratios for 
the combined sample of DLA systems 1 . The Ar/Fe ratio 
is likely to be affected by nucleosynthesis evolution, dust 
depletion and ionization effects. Disentangling these ef- 



1 We follow the usual definition [X/Y] = 
log{N(X l )/N(Y 3 )} - log(X/Y) sun , where i and j indi- 
cate the dominant ionization state in Hi regions for the 
elements X and Y, respectively. In this paper we adopt the 
meteoritic solar abundances of Grevesse & Sauval (1998), with 
the exception of O and Ar, taken from Holweger (2001) and 
SJ98, respectively. 




-2.2 -2.0 -1.8 -1.5 -1.2 -1.0 
[Fe/H] 



Fig. 5. [Ar/Fe] versus [Fe/H] abundance ratios in Damped 
Ly a systems. The sample shown in the figure includes all 
the systems of Table 1 for which a measurement of N (Fe n) 
is available. References to the original measurements are 
given in the table. 



fects is extremely difficult, making the interpretation of 
the data shown in Fig. quite uncertain. As we explain 
below, the ratio of Ar relative to a-capture elements is 
best suited for casting light on ionization effects alone. 



3. The Ar/a ratio in Damped Ly a systems 

In columns 5, 6, and 7 of Table 1 we give the Ar abun- 
dances normalized to those of O, S, and Si, representative 
of a-capture elements measured in DLAs. Most of the sys- 
tems have [Ar/0,Si,S] ratios below the solar value, with 
underabundances of about —0.7 dex in several cases. Only 
the z=3.39 DLA towards Q0000-26 shows solar ratios. The 
[Ar/0,Si,S] ratios show therefore a significant scatter, with 
differences as high as 0.8 dex. All upper and lower limits 
are consistent with this range of values. 

Because O, Si, S, and Ar are all a-capture elements, 
the [Ar/O], [Ar/Si], and [Ar/S] ratios are expected to be 
only weakly dependent on the detailed chemical history 
undergone by each system. In fact, when at least two mea- 
surements of [Ar/O, Si, S] are available for a given absorber, 
one can see in Table 1 that the ratios are in good agree- 
ment between them. The only exception is the z = 3.025 
system in Q0347— 383, for which O and Si abundances 
show a difference of 0.2 dex. However, we do not believe 
that nucleosynthesis is a likely source for explaining the 
significant underabundances found in the [Ar/O, Si, S] ra- 
tios. We now consider both dust and ionization effects as 
possible sources of these underabundances. 



3.1. Dust depletion 
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The evidence that Ar is unlikely to be incorporated in dust 
is mostly based on theoretical arguments (SJ98). In fact, 
local interstellar measurements of A^(Ar i) are feasible only 
in low-column density lines of sight (iV(Hi) < 10 20 cm~ 2 ), 
which are not expected to have a significant dust con- 
tent. The Ar measurements in DLAs provide the possi- 
bility to test the existence of Ar dust depletion in high 
redshift galaxies, in a regime of higher Hi column densi- 
ties. We can estimate the dust content of individual DLAs 
making use of the [Zn/Fe] ratio, given the large differen- 
tial depletion between Zn and Fe, well known from local 
ISM studies (Savage & Sembach 1996) and also supported 
from studies of DLAs (Pettini et al. 1994; Vladilo 1998; 
Prochaska & Wolfe 2002) 2 . In Fig.0we plot [Ar/S] versus 
[Zn/Fe] for the subsample of Table 1 with available Fe and 
Zn data. Should Ar be incorporated into dust, we would 
expect the [Ar/S] ratio to become more and more under- 
abundant with increasing [Zn/Fe], because S is known to 
be undepleted in the interstellar gas (Savage & Sembach 
1996). However, this type of trend is not seen in Fig. 
On the contrary, severe underabundances of Ar are present 
even when [Zn/Fe] is relatively small, i.e. when the amount 
of dust is negligible. This empirical test indicates that a 
mechanism different from dust depletion is responsible for 
the Ar underabundances in DLAs, lending support to the 
theoretical arguments given by SJ98. 

Among the abundances reported in Table 1 only the 
[Ar/Si] ratio is possibly affected by dust depletion, given 
the fact that O and S are known to be undepleted from 
local ISM studies. In the penultimate column of Table 1 we 
give the amount of Si/Ar differential depletion expected 
for DLAs with available Zn data, estimated following the 
procedure of Vladilo (2002). 

3.2. Ionization effects 

Abundance measurements in DLAs are performed using 
dominant ionization states of metals in Hi regions, which 
can be neutral (e.g. N I, O I, and Ar i) or singly ionized 
(e.g. Sill, Sn, Fen, and Znii), depending on the ion- 
ization potential of each element. The predominance of 
such species in Hi regions is well known from local ISM 
studies and can be easily explained in terms of photoion- 
ization equilibrium models of DLAs (Howk & Sembach 
1999; Vladilo et al. 2001, hereafter V&01; Prochaska et 
al. 2002b). Ionization corrections for the abundances can 
be estimated, in the framework of these models, as a func- 
tion of the intensity and spectrum of the adopted ionizing 
continuum. Unfortunately, singly ionized species can also 
arise in intervening Hn regions, if they happen to lie in- 
side the Damped systems. We discuss models with and 
without an intervening Hn region. 
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Fig. 6. [Ar/S] versus [Zn/Fe] abundance ratios in Damped 
Ly a systems. See Table 1 for references to the original 
measurements. 
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2 This is true even if the nucleosynthesis of Zn is not well 
understood (see Prochaska & Wolfe 2002). 



Fig. 7. The Oi/Sin ratio in Damped Ly a systems. 
Diamonds: (from left to right): Levshakov et al. (2002), 
Dessauges-Zavadsky et al. (2001), Molaro et al. (2001). 
Squares: Pettini et al. (2002). Circles: Prochaska et al. 
(2002a). 

3.2.1. Intervening H n regions 

The presence of intervening H n gas in DLAs would yield 
an excess of ionized species, but not of neutrals. As a 
consequence, they would create an (apparent) underabun- 
dance of the Ar/Si and Ar/S ratios, derived by comparing 
An with Sin or Sn. On the other hand, the Ar/O ra- 
tios, derived from neutral species, are not expected to be 
altered by this effect. Therefore, if intervening Hn gas 
exists, we should find [Ar/Si, S] < [Ar/O]. The only two 
[Ar/O] ratios in our sample do not support this possibility 
(see Table 1). Since this sample is too small — and does 
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Fig. 8. The [Ar/S] ratio in Damped Ly a systems. Filled 
diamonds and arrow: measurements and lower limit from 
Table 1; open diamond: z = 2.812 system in Q0528— 2505. 
Continuous curves: predictions of the photoionization 
models S2 (dotted curve) and HI (dashed curve) described 
in Section 3.2.2. 



Fig. 9. The [Ar/Si] ratio in Damped Ly a systems. 
Filled diamonds and arrows: measurements and upper lim- 
its from Table 1; open diamond: z — 2.812 system in 
Q0528— 2505. Continuous curves: predictions of the pho- 
toionization models S2 (dotted curve) and HI (dashed 
curve) described in Section 3.2.2. 



not include cases with significant Ar undcrabundances — 
we have also investigated the behaviour of the Oi/Sin 
ratio, for which 9 DLAs measurements exist, 2 of them 
shown in Table 1 and the others found in the literature 
(Prochaska & Wolfe 1999, Dessauges-Zavadsky et al. 2001, 
Pettini et al. 2002, Prochaska et al. 2002a). The resulting 
[O/Si] ratios are shown in Fig.|7| The mean value of this 
sample is < [O/Si] >= +0.10 ± 0.14 dex. A few Oi mea- 
surements are derived from partially saturated 130.2 nm 
lines, and in these cases the real [O/Si] ratio could be even 
higher. A little excess of [O/Si] may be expected due to dif- 
ferential dust depletion between the two elements (Savage 
& Sembach 1996). Only one system shows a modest un- 
derabundance which can be explained in the framework 
of the ionization balance of H I regions discussed below. 
The general lack of O/Si underabundance indicates that 
intervening Hn gas is uncommon in DLAs. On the ba- 
sis of this result, DLA ionization models which consider 
intervening Hn regions as inherent to DLAs ionization 
structure, such as the models discussed by Izotov et al. 
(2001), are not considered hereafter. 

3.2.2. Ionization balance of the Hi region 

Photoionization models can be constrained when two or 
more ionization states of the same element coexist and 
can be detected in H I regions. This is the case of Al II and 
Aim, the only species of intermediate ionization which 
can be easily measured in DLAs. In fact, Aim shares a 
common velocity distribution — and, hence, a common 
location in the DLAs — with species of lower ionization 



(Wolfe & Prochaska 2000) . General models of DLAs ion- 
ization should be able to explain the abundances of Al ill 
relative to Alii and other ions. In a previous work, a 
trend of decreasing Aim/Alii ratio with increasing N{R i) 
was found and used to constrain photoionization models 
(V&01). Two types of ionizing continuum were consid- 
ered: (1) a hard, QSO-dominated radiation field (Haardt 
& Madau 1996; see also Howk & Sembach 1999) and (2) a 
soft, stellar type radiation field (T e ff — 33, 000 K; Kurucz 
1991). In the first case (HI model) the Aim is produced 
by the hard photons, which are able to penetrate H I layer 
in depth. In the second case (S2 model) the Aim origi- 
nates in a partially ionized interface created by the soft 
radiation field, which is not able to penetrate the neu- 
tral gas in depth. Both types of models are able to repro- 
duce the observed trend Aim/Alii versus A^(Hi) with a 
proper choice of the ionization parameter, U. For model 
HI, the requirement to fit the observed Aim/Alii versus 
A^Hi) trend implies that the ionization parameter must 
vary as U oc A^(Hi)~ 15 . However, the model HI cannot 
be computed at arbitrarily high values of iV(Hi), owing 
to uncertainty in the thermal solution of the photoion- 
ization computations. For model S2 the trend Aim/Alii 
versus TV (Hi) can be reproduced with a constant ioniza- 
tion parameter. We refer to V&01 for more details. After 
constraining the ionization parameter in such a way to 
match the Aim/Alii trend, the HI and S2 models can 
be used to predict the behaviour of At /a ionic ratios as 
a function of iV(Hi). In Figs. [H] and [5] we show the pre- 
dictions for the Ari/Sn and An/Sin ratios, respectively, 
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derived with model HI (dashed line) and S2 (dotted line), 
for a gas with solar abundance ratios. Both models predict 
that deviations from the solar value induced by ionization 
become less important with increasing TV (Hi). The com- 
parison of observations with model predictions in Figs. [SI 
and El yields the following results. 

1. The [Ar/S] and [Ar/Si] ratios, when available for a 
given DLA, give consistent indications on the ionization 
conditions 3 . For instance, both ratios of the z — 3.390 sys- 
tem in QSO 0000-26 are in agreement with model S2. For 
the z = 2A62 system in QSO 0201+365 the ratios yield 
consistent results, slightly above the HI line, when the 
Ar/Si dust correction term is taken into account (Table 
1; dust corrections are negligible in the other cases). The 
concordant indications obtained from different ratios, each 
one measured and modelled independently, indicates that 
the results are, at least, self consistent. 

2. Only a few [Ar/S, Si] ratios lie in the proximity of the 
S2 curve, suggestive of a soft ionizing continuum. Most of 
the data lie instead below the S2 model, with several cases 
close to the HI curve, suggestive of a hard continuum. 

3. Even if the DLA with highest iV(Hi) in our sample 
(the z = 3.390 system in QSO 0000-26) does show negligi- 
ble ionization effects, there is no evidence for the system- 
atic decrease of ionization effects with increasing Af(Hi) 
predicted by the models. The trend could be smeared out 
because different DLAs are embedded in different types of 
ionizing continuum, with one possible case of extremely 
intense continuum, as we mention in the next point. 

4. The z = 2.812 system in QSO 0528-25 is peculiar 
since its redshift is larger than the emission redshift of the 
QSO, z cm = 2.779. Most likely the system has a large, pos- 
itive radial velocity and is very close to the QSO. In this 
case, the large [Ar/a] underabundance at high A^Hi) ob- 
served in this system (open symbol in Figs. |S1 and |5J could 
be due to the effect of a strong, QSO-type radiation field 
enhanced relative to the diffuse background. A detailed 
analysis of this DLA will be presented in a separate paper 
(Peroux et al. 2003). 

4. Redshift evolution of Ar abundances 

We now present evidence for redshift evolution of Ar abun- 
dances, based on measurements in DLAs in the range 
2.3 < z < 3.4 and on local ISM measurements, repre- 
sentative of the conditions at z = 0. 

4.1. Ar abundances in DLAs 

In Fig. 1101 we show all [Ar/a] measurements and limits 
given in Table 1 plotted versus z. For the limits, we adopt 
the most stringent limit on [Ar/a] available. For the mea- 
surements, we define [Ar/a] as the mean value of available 
[Ar/O], [Ar/S], and [Ar/Si] data, with error bars obtained 
from error propagation of individual data. Only for the 

3 The individual DLAs in the figures can be recognized by 
means of their Hi column density, given in Table 1. 
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Fig. 10. The [Ar/a] ratios in DLAs plotted versus redshift 
(bottom axis) and look-back time (top axis), for H = 65 
km s" 1 Mpc -1 , fl m = 0.3, and = 0.7. See Table 1 and 
Section 4. 



z = 3.025 absorber in QSO 0347-38 the individual error 
bars are smaller than the difference between the [Ar/O] 
and [Ar/Si] measurements. In this case we adopt an error 
bar sufficiently large to encompass both measurements. 

The [Ar/a] ratios in Fig. increase with redshift, 
from ~ -0.7 dex at z - 2.4 to ~ dex at z ~ 3.4. The 
trend is supported by all the 10 measurements and the 5 
limits, four of which happen to be stringent. A linear re- 
gression of [Ar/a] versus z for the 10 measurements yields 
a slope m = 0.63 ± 0.12 at 99.92% confidence level (cl). 
If the z = 2.812 system in Q 0528-25 is excluded, the 
linear regression yields m = 0.65 ±0.11 at 99.94% cl. The 
trend versus z is particularly remarkable given the lack of 
any trend versus AT (Hi) as discussed above. This result 
strongly suggests that Ar abundances are subject to evo- 
lution, even if the sample is still small at z > 3 to allow 
firm conclusions to be drawn. In favour of the existence of 
the trend, we note that the linear regression analysis still 
yields consistent results, i.e. m = 0.62 ±0.17 at 98.91% cl, 
after excluding the data point at z = 3.4 (8 measurements; 
z = 2.812 system also excluded). 

To our knowledge, the Ar/a ratio shows the strongest 
redshift variation of relative abundances found so far in 
DLAs. Ionic ratios predicted to be less affected by the 
ionization field do not show such variations. For instance, 
the Si/Fe ratio shows a nearly constant value, although 
with significant scatter, over the large redshift interval 
1.8 < z < 4.2 (Fig. 26 in Prochaska & Wolfe 2002). Once 
corrected for dust effects, the Si/Fe ratio does show some 
increase with redshift (Fig. 3 in Vladilo 2002) which, how- 
ever, can be explained in term of chemical evolution (Fig. 
5 in Calura, Matteucci & Vladilo 2003). As mentioned be- 
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fore, we do not expect significant chemical evolution for 
the Ar/a ratio. The interpretation of the observed evolu- 
tion of Ar/a versus redshift is discussed in Section 5. 

4.2. Interstellar Ar abundances at z = 

Interstellar lines of sight with logV(Hi) > 20.3 through 
the Milky Way or nearby galaxies are representative 
of the behaviour of DLAs at z = 0. Most interstel- 
lar Ar abundances have been obtained for the Milky 
Way. Unfortunately, at the solar metallicity typical of 
the Milky Way, Ar I interstellar absorptions are saturated 
when logTV(Hl) > 20.3 and reliable measurements can be 
obtained only in lines of sight with column densities some- 
what lower than in DLAs. This difference is taken into 
account in the discussion below (Section 5). 

Milky- Way measurements of [Ar/a] are available in 4 
lines of sight, yielding —0.3 < [Ar/O] < 0.0 dex (Jenkins 
et al. 2000; Lehner et al. 2000). Since the local ISM has 
solar metallicity, the Ar underabundance can also be es- 
timated from the [Ar/H] ratios, which are available for 
3 other lines of sight, yielding [Ar/H] = —0.6, —0.4, and 
—0.2 dex (SJ98). Limits obtained from saturated An lines 
are consistent with these values. The median value of all 
available measurements is —0.2 dex, significantly larger 
than the typical [Ar/a] ratio in DLAs at z ~ 2.5. 

The largest Ar deficiencies are found in the Galactic 
plane. Only one, observed towards (3 Cen ([Ar/H] = —0.6 
dex), is comparable to those measured in DLAs at z ~ 2.5. 
According to SJ98, this large deficiency requires the cloud 
to be located close to an ionizing star or to be embedded 
in soft X-ray radiation produced by local hot gas. The ex- 
istence of soft X-ray emitting gas in hot cavities of the 
local H I distribution is well known and the line of sight 
to (3 Cen intersects the Local Bubble (Kuntz & Snowden 
2000; Sfeir et al. 1999) and, probably, a further hot cavity 
associated with Loop I (Iwan 1980; Centurion & Vladilo 
1991). Therefore, the severe Ar deficiency measured to- 
wards (3 Cen may well represent a particular, rather than 
general, interstellar condition. 

Among the local ISM measurements, those taken at 
high Galactic latitude have better chance to reflect the in- 
fluence of radiation fields external to the Milky Way. Only 
three lines of sight have been investigated at \b\ > 35°, 
yielding [Ar/O] = -0.17, -0.18, and -0.27 dex (Jenkins 
et al. 2000; Lehner et al. 2000), consistent with the me- 
dian value of the full sample of Milky- Way measurements, 
-0.2 dex. 

To our knowledge, the only Ar I interstellar absorption 
detected in a nearby galaxy is the one obtained by Vidal- 
Madjar et al. (2000) from a FUSE spectrum of IZwl8 
taken with a large entrance aperture covering the whole 
galaxy. Analysis of the absorption spectrum yields [Ar/Si] 
~ —0.3 dex (Izotov et al. 2001), even though a value some- 
what closer to solar can be derived by taking into account 
the effects of large-scale velocity fields integrated in the 
field of view (Levshakov et al. 2001). 



In summary, with the available data set it seems rea- 
sonable to take [Ar/a] ~ — 0.2 dex as a representative in- 
terstellar value of the local universe. The comparison with 
the low values observed in DLAs at z ~ 2.5 indicates that 
the Ar abundance must significantly increase from z ~ 2.5 
to z = 0. 

5. Discussion 

The present study indicates that the full range of [Ar/a] 
values observed in DLAs can be reproduced in the frame- 
work of photoionization models of H I regions by properly 
tuning the spectrum of the ionizing continuum. Even if 
the models considered are quite idealized, it is clear that 
the magnitude of the underabundance is related to the 
shape of the ionizing spectrum, [Ar/a] < —0.5 dex being 
suggestive of a hard, QSO-type continuum, while [Ar/a] 
closer to solar value being suggestive of a soft, stellar- type 
continuum 4 . The unique capability of An to discriminate 
between different sources of ionization at work in DLAs 
offers the opportunity to cast light on the origin of the ion- 
izing continuum in the associated galaxies, which in prin- 
ciple will include both internal contributions, due to stellar 
emission, and external ones, due to integrated metagalac- 
tic radiation. The origin of the integrated UV background 
has been the object of intensive investigation, with QSOs 
and star-forming galaxies being considered as important 
contributors to the UV background at various redshifts 
(Haardt & Madau, 1996; Songaila 1998; Steidel et al. 2001; 
Bianchi et al. 2001; Haehnelt et al. 2001). From these stud- 
ies a picture is emerging in which the intensity and spec- 
tral shape of the UV background evolves with redshift. 
These studies are necessarily linked to those of the global 
star formation rate in the universe, which also shows evi- 
dence for evolution (Madau, Pozzetti, & Dickinson 1998; 
Lanzetta et al. 2002). Both types of change (UV back- 
ground and star formation rate) are expected to affect 
the ionization state and the Ar abundances in DLAs. 
The trend of Fig. ^| suggests that the ionizing contin- 
uum in DLAs becomes gradually harder from z ~ 3.4 up 
to z ~ 2.3. The interstellar Ar abundances in the local 
universe suggest that the continuum should become softer 
from z ~ 2.3 to z = 0. We discuss the possible implications 
of such evolutionary trends. 

5.1. Ionization properties between z ~ 3.4 and z ~ 2.3 

The general underabundance of Ar at z ~ 2.5 is quite 
clearly established with the present data set. This result 
suggests that the spectrum of the continuum is relatively 
hard, a conclusion which in turn supports an external ori- 
gin of DLAs ionization, with dominant contribution from 
QSOs around redshift z ~ 2.5. This interpretation is con- 

4 Strictly speaking, the latter indication must be taken with 
caution at high column densities since we are not able to com- 
pute the ionization balance for model HI when log iV(H i) > 21 
(see 3.2.2). 
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sistent with the high space density of quasars observed in 
the redshift interval 1.5 < z < 3 (Shaver et al. 1996). 

The interpretation of the modest or negligible under- 
abundances seen at z > 3, suggestive of a softer ionizing 
spectrum, is less straightforward. The decline of space den- 
sity of quasars for z > 3 (Shaver et al. 1996) may lead to 
a softening of the UV continuum if quasars are gradu- 
ally replaced by stars as dominant sources of ionization. 
However, we do not know whether the bulk of the soft 
photons required is internal (e.g. stars inside the DLAs) 
or external (e.g. Lyman-break galaxies). We briefly discuss 
these two possible scenarios, taking also into account the 
possibility that the softening of the ionizing continuum at 
z > 3 may reflect an evolution of the optical depth of the 
IGM. 

UV continuum at z > 3 dominated by internal stars. 
An internal origin of the flux requires a strong synchro- 
nization of stellar emissivity among DLAs, in such a way 
that internal stellar photons dominate at z > 3 and be- 
come negligible at z < 3 in most cases. This requirement 
is rather stringent since it would imply a strong decrease 
of the star formation rate in DLAs at z ~ 3, not easy 
to understand given the weak evidence for redshift evolu- 
tion of metallicities (Pettini et al. 1999, Prochaska & Wolfe 
2000, Vladilo et al. 2000) and relative abundances (Vladilo 
2002) in DLAs. In addition, high redshift DLAs arc likely 
to have low dust content and this makes more critical the 
requirement that the UV continuum is dominated by in- 
ternal starlight. In fact, we know from Milky Way studies 
that the main contribution to the diffuse interstellar radi- 
ation is given by photons scattered by dust (see Bowyer 
2001). If the dust content is low, the UV photons would 
tend to escape from the system rather than ionize the H I 
region(s) inside the system. 

UV continuum at z > 3 dominated by external galax- 
ies. In this scenario the external continuum dominates the 
internal one in the whole redshift interval explored, with a 
spectral distribution gradually changing from hard, QSO- 
type at z < 3 to soft, stellar- type at z > 3. This requires 
(i) the internal emissivity from stars to be generally mod- 
est, (ii) the UV background must become softer at z > 3. 
These requirements seem easy to accommodate within our 
current views of DLAs, LBGs and studies of metagalac- 
tic UV radiation. In fact, the modest stellar emissivity 
of DLAs is consistent with an origin in galaxies charac- 
terized by modest star formation rates, as suggested by 
studies of chemical evolution of DLAs (Calura et al. 2003). 
A low stellar emissivity is also suggested by searches for 
Lya emissions (see Fynbo et al. 1999; Fynbo et al. 2000) 
and H a emissions from DLAs (Teplitz, Malkan & McLean 
1998; Bunker et al. 1999; Bouche et al. 2001) which have 
led, in general, to negative or modest detections. A soft- 
ening of the spectrum of the metagalactic UV continuum 
above z ~ 3 is supported by studies of Siiv/Civ ra- 
tios in metal absorption systems (Savaglio et al. 1997; 
Songaila 1998; see however Kim, Cristiani & D'Odorico 
2002; Levshakov et al. 2002b). If the escape fraction of UV 
photons from Lyman-break galaxies (LBGs) at z > 3 is 



high, as suggested by Steidel et al. (2001), then the change 
of the spectrum could be due to the rise of the contribu- 
tion from LBGs parallel to the decline of the space density 
of quasars. Perhaps more interesting is the possibility that 
the softening of the continuum is related to a change of the 
optical depth of the IGM, such as the one expected dur- 
ing the He n reionization epoch. After this reionization, 
photons with hv > 54 eV are free to propagate through 
the IGM and can easily affect the relative abundance of 
Ar I. In fact, the ratio of the ionization over recombination 
rates is one order of magnitude higher for Ar I than for H I 
at hv > 60 eV (Fig. 3 in SJ98). Evidence for the detec- 
tion of the Hen reionization has been recently reported 
by Bernardi et al. (2002) and Theuns et al. (2002). The 
coincidence between the redshift interval in which the Ar 
abundances change (Fig. I10|) and the interval indicated 
by Theuns et al. (2002) for the He II reionization to occur 
(3.0 < z < 3.4) is indeed quite remarkable. In any case it 
is clear that, from a purely observational point of view, 
the cumulative evidence for a change of the physical con- 
ditions of the IGM around z ~ 3 starts to be compelling. 



5.1.1. Optical depth of the Civ/Siiv layers in DLAs 

The hypothesis that Ar abundances are governed by radi- 
ation fields external to DLAs implies that the H I regions 
must be directly exposed to such radiation. This bears 
implication on the optical depth of the Civ/Siiv layers 
in DLAs, which are kinematically disconnected from the 
layers of low ionization (Wolfe & Prochaska 2000) and 
probably envelope the low ionization gas, in a way similar 
to what observed in the Milky Way, where C iv and Si iv 
layers are stratified at large distances from the Galactic 
disk. In fact, in order to prevent absorption from external 
radiation, the high-ionization layers in DLAs should have 
a negligible Lyman discontinuity. This in turn implies that 
the neutral hydrogen column density associated with the 
highly ionized gas should be iV(Hi) < 10 17 atoms cm -2 . 



5.1.2. Soft-X origin of Ar underabundances in DLAs 

The local ISM studies suggest that strong Ar deficiencies 
might be expected in DLAs if regions of hot, soft X-ray 
emitting gas exist in proximity of the H I region sampled 
by the line of sight. In principle, this possibility could be 
used to interpret the trend seen in Fig. 1101 as entirely due 
to internal sources dominating at all redshifts (stars at 
z > 3 and soft-X bubbles at z <~ 2.5). However, in addition 
to the drawbacks of the "internal scenario" mentioned in 
Section 5.1, this would also require a very large filling 
factor of hot gas in DLAs at 2.4 < z < 3. Since in the 
Milky Way we have detected only one strong Ar deficiency 
of internal origin out of seven measurements, such cases 
do not need to be frequent in DLAs. 
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5.2. Ionization properties between z ~ 2.3 and z = 

The local ISM ionization, representative of the conditions 
at z = 0, is governed by the total contribution of ionizing 
sources internal and external to the Milky Way. The origin 
of the external background in the far UV band is not well 
understood, even though it is clear that (i) some exter- 
nal contribution does exist at high Galactic latitude, (ii) 
sources with soft continuum, such as Galactic starlight 
scattered by dust and integrated light from galaxies, do 
give some contribution; (iii) sources with hard continuum, 
such as integrated QSOs emission and He n Ly a emission 
from the IGM, arc not dominant (Paresce & Jakobsen 
1980; Bowyer 2001; Henry 2002). At low Galactic latitude 
most of the internal background is due to starlight scat- 
tered by dust. Therefore the ionizing continuum in the 
ISM at cither low and high latitudes is expected to be 
soft and, as a consequence, the local Ar deficiencies to be 
small. The modest underabundances of Ar found in the lo- 
cal ISM, with median value —0.2 dex, are consistent with 
this expectation. The fact that local Ar abundances are 
measured in lines of sight with Hi column density lower 
than in DLAs does not affect this result. In fact, these 
lines of sight are less self-shielded and should be more af- 
fected by ionization. A hard continuum would create in 
any case severe Ar underabundances, which are instead 
uncommon. 

If we take the Ar measurements in the local ISM as 
representative of the universe at z = 0, the present results 
indicate that the Ar abundances in DLAs must increase 
from ~ —0.7 dex at z ~ 2.4 to ~ —0.2 dex at z = 0. 
This predicted increase is consistent with the evolution 
of the relative contribution of QSOs and galaxies to the 
UV background in the same redshift interval, derived from 
number density studies of the Lyman a forest. According 
to these studies the contribution of galaxies overcomes 
that of QSOs after z ~ 1 (Bianchi et al. 2001). If this 
is true, we predict to detect a rise of Ar abundances at 
similar, intermediate redshifts. 

6. Conclusions 

From the analysis of the combined sample of Ar abun- 
dance, obtained from An absorptions in DLAs (10 mea- 
surements and 5 limits), we have derived the following 
results. 

Most DLAs show significant Ar underabundances rel- 
ative to other a elements (O, Si, and S), with deficiencies 
as low as ~ —0.8 dex in some cases. 

Departures from cosmic abundances of this magnitude 
are not expected for elements that share a common nucle- 
osynthetic history. We have considered dust and ioniza- 
tion effects as possible explanations of the observed un- 
derabundances. 

From a study of [Ar/a] versus [Zn/Fe] ratios in DLAs 
we have shown that dust depletion does not give a viable 
explanation. This result is consistent with the lack of dust 



depletion expected for Ar from theoretical studies of the 
local interstellar medium. 

As far as ionization is concerned, we have first consid- 
ered the possible effects of intervening H n regions inside 
the DLAs and, on the basis of a study of O i/Si n ratios, we 
have concluded that such H II regions must be uncommon 
and cannot explain the observed underabundances. 

The full range of the observed Ar abundances can in- 
stead be reproduced by photoionization models of H I re- 
gions embedded in an ionizing continuum with variable 
spectral distribution. By using models tuned to match 
the Al m/Al II ratios measured in DLAs wc find that the 
largest deficiencies of Ar can be explained by a hard, QSO- 
dominated continuum, while the modest deficiencies by a 
soft, stellar- type continuum. 

We have found evidence for a redshift dependence of 
the Ar abundances in DLAs. At z ~ 2.5, where most of 
the measurements are concentrated, the deficiencies are 
strong (—0.8/ — 0.6 dex). At z > 3 the deficiencies are 
smaller (> —0.5 dex), with the Ar abundance becoming 
solar at z ~ 3.4. More measurements at z > 3 are required 
to understand how general this latter result is. 

The strong Ar deficiencies at z ~ 2.5 indicate that the 
ionization is dominated by a QSO-type spectrum, which 
we associate with a QSO-dominated metagalactic back- 
ground. The modest deficiencies measured at higher red- 
shifts suggest a predominance of a soft, stellar- type spec- 
trum at z > 3. We have considered an origin of this soft 
spectrum in both starlight internal to DLAs and galac- 
tic emission external to DLAs. In both cases, the gradual 
change of the ionizing continuum from z ~ 2.4 to z ~ 3.4 
poses strong requirements on the nature of DLAs and/or 
the origin of the metagalactic background. 

If the redshift variation of Ar abundances is due to 
the evolution of internal, stellar emission, it requires a 
synchronization of evolution of DLAs, with most systems 
having strong star formation rates at z > 3 and weak at 
z < 3. This seems to be at odd with the lack of clear signal 
of evolution of DLAs abundances. Also the need for dust 
at very high redshift, required to maintain an internal, 
diffuse source of UV radiation, seems to be difficult to 
reconcile with the low dust content expected in the early 
evolutionary stages of DLAs. 

If the trend is induced by the evolution of the meta- 
galactic ionizing background, it requires (i) the H I regions 
in DLAs to be directly exposed to the external back- 
ground, (ii) the contribution from internal emission to 
be modest, and (iii) the external background to become 
softer at z > 3. The first requirement would imply that 
the Siiv/Civ layers associated to DLAs should be opti- 
cally thin to ionizing radiation. The second requirement 
is consistent with previous claims of low star-formation 
rates based on abundance studies and searches for Lya 
and H a emission in DLAs. The third requirement is con- 
sistent with the evolution of the metagalactic spectrum 
inferred from Siiv/Civ measurements in the IGM, with 
claims that LBGs at z > 3 may have a high escape fraction 
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of ionizing photons, and with the recent finding that the 
He II reionization does occur between z ~ 3.4 and z ~ 3.0. 

We have summarized the results of studies of Ar abun- 
dances in the local ISM taken as representative of a DLA 
system at z = 0. The typical underabundance of Ar is 
~ —0.2 dex, consistent with the expectation that the ion- 
ization in the local universe is dominated by sources with 
soft continuum. Only in one case a relatively strong un- 
derabundance of Ar is found, which may originate in gas 
embedded in hot, soft-X ray emitting gas known to ex- 
ist in the local ISM. Taken together, the measurements of 
Ar abundances in DLAs and in the local universe appear 
to offer a new tool for probing the redshift evolution of 
the ionization conditions in the universe. Present obser- 
vational limitations prevent us to use this tool to probe 
intermediate redshifts for determining the epoch at which 
QSOs stop being dominant contributors to the metagalac- 
tic UV background. This epoch is estimated to be at z ~ 1 
from recent studies of the Ly a forest. An independent as- 
sessment based on Ar abundances will be possible from 
observations with UV spectrographs fed by space-born 
telescopes with large collecting areas filling the gap of Ar I 
measurements between z = and z ~ 2.3. 
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